Transverse Spin-Orbit Force in the Spin Hall Effect in Ballistic Semiconductor Wires 
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We introduce the spin and momentum dependent force operator which is defined by the Hamil- 
tonian of a clean semiconductor quantum wire with homogeneous Rashba spin-orbit (SO) coupling 
attached to two ideal (i.e., free of spin and charge interactions) leads. Its expectation value in 
the spin-polarized electronic wave packet injected through the leads explains why the center of the 
packet gets deflected in the transverse direction. Moreover, the corresponding spin density will be 
dragged along the transverse direction to generate an out-of-plane spin accumulation of opposite 
signs on the lateral edges of the wire, as expected in the phenomenology of the spin Hall effect, when 
spin-t and spin- j polarized packets (mimicking the injection of conventional unpolarized charge cur- 
rent) propagate simultaneously through the wire. We also demonstrate that spin coherence of the 
injected spin-polarized wave packet will gradually diminish (thereby diminishing the "force") along 
the SO coupled wire due to the entanglement of spin and orbital degrees of freedom of a single 
electron, even in the absence of any impurity scattering. 

PACS numbers: 72.25.Dc, 71.70.Ej, 03.65.Sq 



The classical Hall effect 1 is one of the most widely 
known phenomena of condensed matter physics because 
it represents manifestation of the fundamental concepts 
of classical electrodynamics — such as the Lorentz force — 
in a complicated solid state environment. A perpendicu- 
lar magnetic field B exerts the Lorentz force F = qv x B 
on current I flowing longitudinally through metallic or 
semiconductor wire, thereby separating charges in the 
transverse direction. The charges then accumulate on 
the lateral edges of the wire to produce a transverse "Hall 
voltage" in the direction gl x B. Thus, Hall-effect mea- 
surements reveal the nature of the current carriers. 

Recent optical detection^ of the accumulation of spin- 
t and spin- 1 electrons on the opposite lateral edges of 
current carrying semiconductor wires opens new realm of 
the spin Hall effect This phenomenon occurs in the ab- 
sence of any external magnetic fields. Instead, it requires 
the presence of SO couplings, which are tiny relativistic 
corrections that can, nevertheless, be much stronger in 
semiconductors than in vacuum^ Besides deepening our 
fundamental understanding of the role of SO couplings 
in solids, 4,5 the spin Hall effect offers new opportunities 
in the design of all-electrical semiconductor spintronic 
devices that do not require ferromagnetic elements or 
cumber some-to-control external magnetic fields. 5 

While experimental detection of the strong signatures 
of the spin Hall effect brings to an end decades of theoret- 
ical speculation for its existence, it is still unclear what 
spin-dependent forces are responsible for the observed 
spin separation in different semiconductor systems. One 
potential mechanism — asymmetric scattering of spin-| 
and spin- 1 electrons off impurities with SO interaction — 
was invoked in the 1970s to predict the emergence of 
pure (i.e., not accompanied by charge transport) spin 
current, in the transverse direction to the flow of longi- 
tudinal unpolarized charge current, which would deposit 
spins of opposite signs on the two lateral edges of the 
sample. 6 However, it has been argued 7 that in systems 
with weak SO coupling and, therefore, no spin-splitting 



of the energy bands such spin Hall effect of the extrinsic 
type (which vanishes in the absence of impurities) is too 
small to be observed in present experiments 2 (unless it 
is enhanced by particular mechanisms involving intrinsic 
SO coupling of the bulk crystal^). 

Much of the recent revival of interest in the spin Hall 
effect has been ignited by the predictions^^ for substan- 
tially larger transverse pure spin Hall current as a re- 
sponse to the longitudinal electric field in semiconductors 
with strong SO coupling which spin-splits energy bands 
and induces Berry phase correction to the group velocity 
of Bloch wave packets. 11 However, unusual properties of 
such intrinsic spin Hall current in infinite homogeneous 
systems, which depends on the whole Fermi sea (i.e., it 
is determined solely by the equilibrium Fermi-Dirac dis- 
tribution function and spin-split Bloch band structure) 
and it is not conserved in the bulk due to the presence 
of SO couplingj2ii£ have led to arguments that its non- 
zero value does not correspond to any real transport of 
spins 12,13 so that no spin accumulation near the bound- 
aries and interfaces could be induced by any intrinsic 
mechanism (i.e., in the absence of impurities^ 3 -) . 

On the other hand, quantum transport analysis of spin- 
charge spatial propagation through clean semiconductor 
wires, which is formulated in terms of genuine nonequi- 
librium and Fermi surface quantities (i.e., conserved spin 
current giii^i^ and spin densitiesii) , predicts that spin 
Hall accumulation 2,3 of opposite signs on its lateral edges 
will emerge due to strong SO coupling within the wire re- 
gion. 17 Such mesoscopic spin Hall effect is determined by 
the processes on the mesoscale set by the spin precession 
length^i 7 - and depends on the whole measuring geom- 
etry (i.e., boundaries, interfaces, and the attached elec- 
trodes) due to the effects of confinement on the dynamics 
of transported spin in the presence of SO couplings in 
finite-size semiconductor structures 

Thus, to resolve the discrepancy between different the- 
oretical answers to such fundamental question as — Are 
SO interaction terms in the effective Hamiltonian of a 
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clean spin-split semiconductor wire capable of generating 
the spin Hall like accumulation on its edges? — it is highly 
desirable to develop a picture of the transverse motion of 
spin density that would be as transparent as the famil- 
iar picture of the transverse drift of charges due to the 
Lorentz force in the classical Hall effect. Here we offer 
such a picture by analyzing the spin- dependent "force", 
which can be associated with any SO coupled quantum 
Hamiltonian, and its effect on the semiclassical dynamics 
of spin density of individual electrons that are injected as 
spin-polarized wave packets into the Rashba SO coupled 
clean semiconductor quantum wire attached to two ideal 
(i.e., interaction and disorder free) leads. 

The effective mass Hamiltonian of the ballistic Rashba 
quantum wire is given by 

* 2 

# = ^ + ^(<xxp)-z + Vconiiy), (1) 

where p is the momentum operator in 2D space, & = 
( y & x , , <r z ) is the vector of the Pauli spin operators, and 
Konf (y) is the transverse potential confining electrons to 
a wire of finite width. We assume that the wire of dimen- 
sions L x x L y is realized using the two-dimensional elec- 
tron gas (2DEG), with z being the unit vector orthogonal 
to its plane. Within the 2DEG, carriers are subjected to 
the Rashba SO coupling of strength a, which arises due 
to the structure inversion asymmetry 4 (of the confining 
potential and differing band discontinuities at the het- 
erostructure quantum well interface 2 ^) . 

This Hamiltonian generates a spin-dependent force op- 
erator which can be extractedSii 2 ^ within the Heisenberg 
picture^ 3 - as 

Fh = m *^r = W [HA * H ' H]] (2) 

2<* 2 m* dV cont (y H ) 

= — (p^^%- dilH y- 

Here the Heisenberg picture operators carry the time 
dependence of quantum evolution, i.e., Pij(t) = 

e iHt/h£ e -iHt/h^ & ^ = e iHt/h & z e -iHt/n^ and = 

e iHt/hjj e -iHt/h ^ w h ere g-z ^ p and - are ^ n ^ e Schrodinger 
picture and, therefore, time-independent. 

Since the force operator 22 depends on spin through 
a z H , which is a genuine (internal) quantum degree of free- 
dom, 23 it does not have any classical analog. Its phys- 
ical meaning (i.e., measurable predictions) is contained 
in the quantum- mechanical expectation values, such as 
(F y )(t) = (®(t = 0)\F%(t)\y(t = 0)) obtained by acting 
with the transverse component Fjj of the vector of the 
force operator (F§, Fjj) on the quantum state \ty(t = 0)) 
of an electron. While such "force" can always be asso- 
ciated with a given quantum Hamiltonian, its usefulness 
in understanding the evolution of quantum systems is 
limited — the local nature of the force equation cannot 
be reconciled with inherent non-locality of quantum me- 
chanics. For example, if the force "pushes" the volume of 



a wave function locally, one has to find a new global wave 
function in accord with the boundary conditions at infin- 
ity (the same problem remains well-hidden in the Heisen- 
berg picture where time dependence is carried by the op- 
erators while wave functions are time-independent). Nev- 
ertheless, analyzing the dynamics of spin and probability 
densities in terms of the action of local forces can be in- 
sightful for particles described by wave packets (whose 
probability distribution is small compared to the typical 
length scale over which the force varies)™* 2 ^ 

Therefore, we examine in Fig. the transverse SO 
"force" (F y ) in the spin wave packet state, which at t = 
resides in the left lead as fully spin-polarized (along the 
z-axis) and spatially localized wave function 2 ^ 2 ^ 

9(t = 0) = Csin ( — ny ) e ik * x ~ 6k l x2 /* Xa . (3) 
\{Ly + l)aj 

This is a pure and separable quantum state = 0)) = 
|$) <g> \cr) in the tensor product of the orbital and spin 
Hilbert spaces Ti 7Y S . The orbital factor state (x, y\<&) 
consists of the lowest subband of the hard wall transverse 
confining potential and a Gaussian wave packet along the 
x-axis whose parameters are chosen to be k x a = 0.44 and 
Sk x a = 0.1 (C is the normalization constant determined 
from ($|$) = 1). The spin factor state is an eigenstate 

of g z , i.e., xr = ^ q ^ o r Xi = ^ \ ^ • 

Unlike the case 21 of an infinite 2DEG, the exact so- 
lutions of the Heisenberg equation of motion for c> z H {t\ 
yii(t) and p#(£) entering in Eq. (0 are not available 
for quantum wires of finite width. Thus, we compute 
the expectation value (Sfr(t)\F y \i&(t)) in the Schrodinger 
picture by applying the evolution operators e - %Ht l h 
present in Eq. (0) on the wave functions \^(t)} = 
Y.n e ~ iEnt/h \ E n)(F n \^(t = 0)). To obtain the exact 
eigenstate a 24 i 26 i 27 \E n ) and eigenvalues E ni we employ 
the discretized version of the Hamiltonian Eq. Q}. That 
is, we represent the Hamiltonian of the Rashba spin-split 
quantum wire in the basis of states |m) (g> |cr), where 
|m) are s-orbitals (r|m) = ip(r — m) located at sites 
m = (Tn x ,m y ) of the L x x L y lattice with the lattice 
spacing a (typically 19 a ~ 3 nm). This representation 
extracts the two energy scales from the Rashba Hamil- 
tonian Eq. QJ: t Q = h 2 /(2m*a 2 ) characterizing hopping 
between the nearest-neighbor sites without spin- flip; and 
tso = ct/2a for the same hopping process when it involves 
spin flip.iiiS The wave vector of the Gaussian packet 
k x a = 0.44 is chosen 17 to correspond to the Fermi energy 
Ep = — 3.8t Q close to the bottom of the band where tight- 
binding dispersion relation reduces to the parabolic one 
of the Hamiltonian Eq. QJ. In this representation one 
can directly compute the commutators in the definition 
of the force operator Eq. thereby bypassing subtleties 
which arise when evaluating the transverse component of 
the force operator —dV con f(yH)y/dyH stemming from the 
hard wall boundary conditions. 28 

Figure Q] shows that as soon as the front of the spin- 
polarized wave packet enters the strongly SO coupled re- 
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FIG. 1: (Color online) The expectation value of the trans- 
verse component of the SO force operator (upper panel) in 
the quantum state of propagating spin wave packet along 
the two-probe nanowire. The middle panel shows the cor- 
responding transverse position of the center of the wave 
packet as a function of its longitudinal coordinate. The ini- 
tial state in the left lead is fully spin-polarized wave packet 
Eq. which is injected into the SO region of the size 

L x x L y = 100a x 31a (a ~ 3 nm) with strong Rashba cou- 
pling tso = a /2a = 0.1 t and the corresponding spin preces- 
sion length Lso = 7r£ a/2£so ~ 15.7a < L x (middle panel) or 
weak SO coupling tso = 0.01£ o and Lso ~ 157a > L x (lower 
panel) . 



gion, its center (y)(t) = (® (t)\y\S& (t)) will be deflected 
along the ?/-axis in the same direction as is the direc- 
tion of the transverse SO "force". However, due to its 
inertia the packet does not follow fast oscillations of the 
SO "force" occurring on the scale of the spin precession 
lengthiLiS Lso = 7rt a/2tso on which spin precesses by 
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FIG. 2: (Color online) The dynamics of spin density S(r) = 
[S x (r), S y (r), S z (r)] induced by simultaneous propagation 
of two electrons through quantum wire 100a x 31a with the 
Rashba SO coupling t so = 0.1t o . Both electrons are injected 
at t — from the left lead, one as spin-| and the other one as 
spin-j polarized (along the z-axis) wave packet Eq. @. The 
different snapshots of the sum of their spin densities are taken 
at the points (a), (b), (c) where the transverse SO "force" and 
the y- coordinate of the center of these wave packets have val- 
ues shown in the middle panel of Fig. Q 



an angle tt (note that the spin splitting generates a finite 
difference of the Fermi momenta, which is the same for 
all subbands of the quantum wire in the case of parabolic 
energy-momentum dispersion, so that Lso is equal for all 
channels 26 ). In contrast to an infinite 2DEG of the intrin- 
sic spin Hall effect in quantum wires electron mo- 
tion is confined in the transverse direction and the effec- 
tive momentum-dependent Rashba magnetic field B#(k) 
is, therefore, nearly parallel to this direction™^ Thus, 
the change of the direction of the transverse SO "force" 
is due to the fact that the z-axis polarized spin will start 
precessing within the SO region since it is not an eigen- 
state of the Zeeman term a • B#(k) [i.e., of the Rashba 
term in Eq. ([T])]. 

The transverse SO "force" and the motion of the cen- 
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ter of the wave packets in Fig. [I] suggests that when two 
electrons with opposite spin-polarizations are injected si- 
multaneously into the SO coupled quantum wire with 
perfectly homogeneous 25 Rashba coupling, the initially 
unpolarized mixed spin state will evolve during propaga- 
tion through the wire to develop a non-zero spin density 
at its lateral edges. This intuitive picture is confirmed 
by plotting in Fig. the spin density, 



S m (t) = -<*(t)|<r®|m><m|*(t)> 



(4) 



corresponding to the coherent evolution of two spin wave 
packets, |tf(t = 0)) = |*) (8) |t> and \9(t = 0)) = |$)<8>||), 
across the wire. 

The mechanism underlying the decay of the transverse 
SO "force" intensity is explained in Fig. |3 where we 
demonstrate that (initially coherent) spin precession is 
also accompanied by spin decoherencei^^ These two 
processes are encoded in the rotation of the spin polariza- 
tion vector P and the reduction of its magnitude (|P| = 1 
for fully coherent pure states p 2 s = p 3 ), respectively. The 
spin polarization vector is extracted from the density ma- 
trix p s = (1 + P • <x)/2 of the spin subsystem. 23 The spin 
density matrix p s is obtained as the exact reduced den- 
sity matrix at each instant of time by tracing the pure 
state density matrix p(t) = \i& (t)) (Sfr (t)\ over the orbital 
degrees of freedom, 

p s {t) = Tr \*(t)){*(t)\ = J2(™\*(t))(*(t)\m) 

m 

= E c m , a {t)\a){a'\c* mt At)- (5) 
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FIG. 3: (Color online) Spin precession, as signified by oscil- 
lations of the spin polarization vector (P x , P y , P z ), and spin 
decoherence (as measured by decrease of the purity |P| below 
one) of the spin state of a single electron propagating along 
the Rashba quantum wire 100a x 31a with the SO coupling 
strength tso = 0.1t o (Lso ~ 15.7a). The electron is injected 
from the left lead as a spin-| polarized wave packet, whose 
spin subsystem is therefore fully coherent |P| = 1 at t = 0. 
The bottom panel shows the z-component of the spin density 
S z (r) at different values of \P Z — J dr S z (r) (selected in the 
upper panel) along the wire. 



The dynamics of the spin polarization vector and the spin 
density shown in Fig. 01 are in one-to-one correspondence 



-P(t) = -Tr s =£)S m (t). (6) 

The incoming quantum state from the left lead in Fig. 
is separable = 0)) = T, m ,* c ™,*(t = 0)|m) \a) = 
|$) <S> | T)j an ^ therefore fully spin coherent |P| = 1. 
However, in the course of propagation through SO cou- 
pled quantum wires it will coherently evolve into a non- 
separable 23 state where spin and orbital subsystems of 
the same electron appear to be entangledw^SL Note that 
Fig. also shows that at the instant when the center of 
the wave packet enters the wire region, its quantum state 
is already highly entangled as quantified by the non-zero 
von Neumann entropy (associated with the reduced den- 
sity matrix of either the spin p s or the orbital subsystem 



Po) 



S{p s ) = S{p ) = - 




which is a unique measure^ of the degree of entangle- 
ment for pure bipartite states (such as the full state \^(t)) 
which remains pure due to the absence of inelastic pro- 
cesses along the quantum wire). 

While this loss of spin coherence (or polarization) is 
analogous to the well-known DP spin relaxation in diffu- 
sive SO coupled systems here the decay of the spin 
polarization vector takes place without any scattering 
off impurities (or averaging over an ensemble of elec- 
trons propagating through ballistic SO coupled quantum 
dot structures 18 ). Instead, it arises due to wave packet 
spreading (cf. lower panel of Fig.[3|), as well as due to the 
presence of interfaces 19 (the wave packet is partially re- 



5 




FIG. 4: (Color online) The spin accumulation (S x (r), S y (r), 
S z (r)) induced by the ballistic flow of unpolarized charge cur- 
rent, simulated by injecting one after another 600 pairs of 
spin- 1 and spin-j polarized (along the z-axis) wave packets 
from the left lead, through quantum wire 100a x 31a with the 
Rashba SO couplings: (a) tso = 0.1t o (Lso ~ 15.7a) and (b) 
tso = 0.01t o (Lso « 157a). 



fleeted at the lead/SO-region interface for strong Rashba 
coupling) and boundarie s 18 ! 19 of the confined structure. 
Thus, the decoherence mechanism revealed by Fig. |3| is 
also highly relevant for the interpretation of experiments 
on the transport of spin coherence in high-mobility semi- 
conductor^ and molecular spintronic devices. 33 

The interplay of the oscillating and decaying (induced 
by spin precession and spin decoherence, respectively) 
transverse SO "force" and wave packet inertia leads to 
spin-j electron exiting the wire with its center deflected 
toward the left lateral edge and the spin-j density appear- 
ing on the right edge 17 for strong SO coupling tso = 0.1t o 
in Figs, n an d [21 This picture is only apparently coun- 
terintuitive to the naive conclusion drawn from the form 
of the force operator itself Eq. (0), which would suggest 
that spin- 1 electron is always deflected to the right while 
moving along the Rashba SO region. While such situ- 
ation appears in wires shorter than Lso (as shown in 
the lower panel of Fig. [Q, in general, one has to take 
into account the ratio L x /Lso, as well as the strength 
of the SO force oc a 2 , to decipher the sign of the spin 
accumulation on the lateral edges and the sign of the 
corresponding spin currents that will be pushed into the 
transverse leads attached at those edges 

When we inject pairs of spin-j and spin-j polarized 
wave packets one after another, thereby simulating the 
flow of unpolarized ballistic current through the lead- 
wire-lead structure (where electron does not feel any 
electric field within the clean quantum wire region) ^£ we 
find in Fig. 0] that the deflection of the spin densities of 



individual electrons in the transverse direction will gen- 
erate non-zero spin accumulation components S z (r) and 
S x (r) of the opposite sign on the lateral edges of the wire. 
While recent experiments find S z (r) with such properties 
to be the strong signature of the spin . Hall effect, 2,3 here 
we confirm the conjecture of Ref. H7I that S x (r) can also 
emerge as a distinctive feature of the mesoscopic spin Hall 
effect in confined Rashba spin-split structures — it arises 
due to the precession (Fig. |3J) of transversally deflected 
spins. Note that S x (r) ^ accumulations cannot be 
explained by arguments based on the texturelike struc- 
ture 2 ^ of the spin density of the eigenstates in infinite 
Rashba quantum wires wher o 26 i 27 S x (r) = 0. 

In conclusion, the spin-dependent force operator, de- 
fined by the SO coupling terms of the Hamiltonian of 
a ballistic spin-split semiconductor quantum wire, will 
act on the injected spin-polarized wave packets to deflect 
spin-j and spin-j electrons in the opposite transverse di- 
rections. This effect, combined with precession and deco- 
herence of the deflected spin, will lead to non-zero z- and 
x-components of the spin density with opposite signs on 
the lateral edges of the wire, which represents an example 
of the spin Hall effect phenomenology^*^ that has been 
observed in recent experiments^^ The intuitively ap- 
pealing picture of the transverse SO quantum- mechanical 
force operator (as a counterpart of the classical Lorentz 
force), which depends on spin through cP, the strength of 
the Rashba SO coupling through a 2 , and the momentum 
operator through the cross product p x z, allows one to 
differentiate symmetry properties of the two spin Hall ac- 
cumulation components upon changing the Rashba elec- 
tric field (i.e., the sign of a) or the direction of the packet 
propagation: S z (r)_ a = S z (r) a and S z (r)_ p = -S*(r) p 
vs. S x (r)- a = —S x (r) a (due to opposite spin preces- 
sion for —a) and S x (r)- P = S x (r) p . These features are 
in full accord with experimentally observed behavior of 
the S z (r) spin Hall accumulation under the inversion of 
the bias voltage, 3 as well as with the formal quantitative 
quantum transport analysis 17 of the nonequilibrium spin 
accumulation induced by the flow of unpolarized charge 
current through ballistic SO coupled two-probe nanos- 
tructures. 

Finally, we note that a 2 dependence of the transverse 
SO "force" is incompatible with the a-independent (i.e., 
"universal") intrinsic spin Hall conductivity <j s h = c/Stt 
(describing the pure transverse spin Hall current j* = 
(TshEx of the z-axis polarized spin in response to the 
longitudinally applied electric field E x ) of an infinite ho- 
mogeneous Rashba spin-split 2DEG in the clean limit, 
which has been obtained within various bulk transport 
approaches JSiiiiiSii^i On the other hand, it supports 
the picture of the SO coupling dependent spin Hall accu- 
mulations- 7 - S z (r), S x (r) and the corresponding spin Hall 
conductances^ (describing the z— and the x-component 
of the nonequilibrium spin Hall current in the trans- 
verse leads attached at the lateral edges of the Rashba 
wire) of the mesoscopic spin Hall effect in confined struc- 
turespi^^ By the same token, the sign of the spin ac- 
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cumulation on the edges (i.e., whether the spin current 
flows to the right or to the left in the transverse direc- 
tion 15 ) cannot be determined from the properties 34 of 
<j s h . Instead one has to take into account the strength of 
the SO coupling a and the size of the device in the units 
of the characteristic mesoscale Lso, as demonstrated by 
Figs, n and 21 This requirement stems from the os- 
cillatory character of the transverse SO "force" brought 



about by the spin precession of the deflected spins in the 
effective magnetic field of the Rashba SO coupled wires 
of finite width. 
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